


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1993-03 


Resolution in radar mapping 


Anderson, Michael D. 


Monterey, California: Naval Postgraduate School 
http://ndl.handle.net/10945/24223 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 


f (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist : Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

| | LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


. heey an a 
* = 
C a re 


rie) 


a Ld 
er 4 PORTE SLi 
* 


he led oe Cd 
te €pP amet 2.9 
peal see TTT ha 


taut inate 
A a a8 c 
aed eee et lk . 7 3 ; ' 
e & Spree * e > = - * i J : iy ‘ 
a | : 5 i 4 ul a e 7 
\- % LD) Py aed dace 
ye Petey 
ar) ‘ Pe Per bs ti 3 
a t . : ri asa Pao 
Fy i % 
‘ ’ 
CT tid el ind thd M4 ar oy, | 
Be o Ree ete ee teed 
ee! 


o oem 8G Bade ry - ’ 
q ; 
La ls et iy y 
‘ ; rl ss) 
oF DL oh Bi aided 
Par ie ete 
cy 


Pt 
i i } ey ae SCRAPE re ett 
+ are roe Aer Pa 
ae! 


ws taf. 
araddeslate ytat 
Oa dates “ Freee 
ae Py or 4 oy 
y 


ar * % Hy 
Tad A 
re] t f he 
3 Be, Oo bp at 
i Cras 


an 
vr 
Sonu F 
CRETE ra 7 LT 4 ab 
SEN : F ie o oe tababste Se 
Ane A ie hee Ry ee ie or 
a? p , eS | pet 
% Pr ad ed 
pepe) 
. tra 
Ca inl 


D 
Orr a ees 
ote apts tty x A F : 
a4 7 1 5 -) a 
o. Die af Fe > eo? e A | er H 
Pepe a oc Sh vey ok 
rans ins 1] bist i , *. ae : 
a ee) t . be er 
j : Larter try Ti 
ai 


Phd 

i ast gt tetas 
Py he oe A : “ TH ae re 

Sar i i UE bed Pi b 2 f 7 PTreti Te iy Ae H f “Ot +] 

vat hr ‘ Ma e«aSf 
: x 4 : ° : rar iF maT Por Hand adh vg : 
eh A 0 , on 3 
‘ 3 By tear biers br Oak 
are? s s 

ae 

rs 


F 
Pray 


Pel CY o fs . 3 
: tT 
H h 
eS) Th 
Paty 
y 


rat ee 
gal * > 
© pee tte rd 
rs ee ee) te 0d 7 

Po LY OL ee p ‘ , A ee 
ry i aM ay % rc Pere Prva s ba eke 

: : F) ; ' A res cas ee 

Pra) 2 auytad 


ebeu Vike tr 
pee age tothe bsbs 95 r : 3 
Pi Fi ory BL ORH Pant. | Oe Or , c M 
5 ; Pa thea Peel 4 L a { J 
Pritts SP in Ser YY Ae ae / PH iad EP) pra edleate ; wlateegse: 
rie mea eee Loe Sa fel bt Cre? Ae hed - | rvs ee a Br 
a arr Be Paar x Pa ae ‘ae per ee ely teed 
rad ay A 4 rept Fee y P rrr) a Pee 

: PEER Mee er aa, Sia) ; Cer ¢ ca eet hee by Sak oe poe SF aE ie Bee me eer 

rw ; a ar Pras oT A ie : P ria ac E : ; aT ds 4 J = ota tgtal aS lar id ae te Or Edt 
d al 2 ’ ae : “ a, "y k 5 ra Oh be Peay th hid Be OF Soeey 
* . rl Pat rs Z 4 + Py Pd 
ts ; URI H NA AE ; cape hah 
j ee J ea b ) L Cs . Dect Dip So 

A 

- 7 ous > a ra ee? 


sj 
‘ Pe Lyi! 


“era hale ® 
Pia a iF 
. A 4 
a fy TY ‘> k Te 3 b = LE 
tnd Sov 5 ah ae 5 
Prod * rR . ad be P A a 
TT * iar et an PS | ; en F ’ Se HPD, 
ri a aS i oie a “ A . 4 yy i ee. Tah PATEL | wa 
oe Ps ae ; Maid Poftsdadt th 
F Park Lia it ees “ 
f 


— 


py panes Any eT 
P-t r e eagere E 
m, xt, 
. 


Prey 
ie PL ta 


BP ettree 8 


a 
Pr re a eee 
pe Oe Tae Zz 
wee. 


LW hhh 
Pi Paes Set bt 4 
RTT tr eo ft, 
oO 4r™ ee Dt deel tae 
Perri 
Pe Fi Beh y 


Me rag 
o 4° a4 


G 


oe ‘ 
re LT hl ad J 
Set te agentes 
Peet atk Meus 
hapa oy Ifa) ae E 
Pay 
bat Pe y 
afd mae yi 
c a a ae 
aah Mila tee Bi 


ad 
i 


¥ ] bead aL 
PCT, ial ee a ea stve rae hed eee 
OE ne ay ie ate 


i 
Pee es 


t- 
Pane , Sra 
meg : . ri . F ; ae 
yee c Pa tae 
A Cater 
yttaytoect tae 
Sy 
. { fi 4 Pee. 


ar Leg 
3 : rot tae 
i Tol od 


e 


’ 
Phi 
etgtnee 


' 
, are 

my & A Pes eren yr 

E wh get args 

Cater fete 


oy 


“a5 
eer a 


ae 


at 
Le] 
ero 


Or en ee ee 


met "att ee 
Pee Ee ty Te 
gd 


. 
~ 


as 


sadeneni a 
rr ar er a t 3 
Pee ray | . A 
A A : 
i aie Pre 
bole nie 
Te i | Pel dal et aL fet 
FARCE Oe Pr ga ah Pha ergy y 
« OUT a Be tr os ree a 
“aye? r) 7 rie by: 
rrr ore ks 
oe Maren yg yee ae 
Pee ae 
rd 
eee anree! 
ye ee Pe f me Pe 
. 1823) peer e 
aici asters 
reste 


Ch tag 1 4a Hi A by cee 4 
CH ot oi oy ih Ae tae I ae ny a sy re : 
x 5 (s iene me yr es Sige . Se adeeebe et te eo 
Vy ee Bar) 4 Y Fahad hie Pr eRe TS 
Rr La eee ¥ t realy 2a f Po ad ie ote 
Po 

hi ed = 4 zg 
. , arg ye” A . 4 pe 
pecreere te ‘ ee Sor ES AP ae ? 
+ dab hd teh EE stepped 
sgetces an Atte F ee nae 
es p # Sr et et cate feo 
Ki ee Tiel ee ah nae es 
ee tte tel oe 


arto) 
caeeaniws 
Fer a 
ef 


eer see . H 7 
ater ees sorte a 4 5 mf 7 ; ae ‘ me a Fe 
a Ora abe rad pee ea 


P 
eee ii st a he Kad, 
thas parr , het L hl bid at 25 fe 
we == ent i Ca§ 

CLs quoter tae 
Potties 

rk Loe | rar ad 

‘3 


nD preted 5 
re ee Ld rr os 
weer tese ; Picea pees 2 oar tae 
ares ERTS one rtae ame PN 
Qs Ba Yea aa ae > tes te = 
Pry Deter rr TE POT CTL hee cae et vat ssa : ri A ‘ papa trae a 
Pare y pes 7 
in HY ea lf Hi est i r $i eats a eres Pa dideed 7 Pt enti bY) 1 Cl deb ihaw ali eee nae] =f 
eater 4 arts ] ‘ eR hd LIND Le er ee 5 4 pair ee We ad eta del 
Ce Pa ty Mey Ae ate 2 ig tty aie id r eats aa GA Ce Rig inne ng 4 
ee, gt athe ey ‘ie . eae ee os ed ve ity sr satis] : ee 
Oe TD r , POrPreT ree LAA 2 “4 £: 
is arab ‘ FA rr array oy ergo a Prd Re ‘ a ra betel eet 
TP ir 5 Thales) : 5 j gtatet ji nt Pte ee ete 
oye pee iierioeh Et eet Has errata Seba Brie) i eh seestbts ais 
Senet eae a PROT Tester Pt 
; ‘ a Ot aed 4 Pen he ee os ; a 
Ae : ; 0 tit 3 ube 
a hale ae 1%y ayy + eerie ty Pre i fierce 
oY Fi a _ tees Oy ba Ae eenaece Fl 
re c i ot’ Hg f 5 Ff Brrr a be hey ’ f 1 
Hf Soe ae *y ‘ * aa ris ea et 
SEU batt ! vie ay : A $ i rere La CEP tiie 
OR Ree Aa Rat Bet e TS 
aye eta! ee rate che ae ee 


so 
a 


pee BESS 


* Tpen? pra tpt. *) ee ee! Fs 

as ! Sait z 7 Ce . 

A n > ot yrat i ; eae rt} alas 

ri eet ha be 2 BS ‘hia , ort a 3 mi creat 01 FS . ae ACAD Te stiri aH es 

Pore Mer 5 ai pri ite. 7 : fi La] pes iJ stared * eens eT te 

A } Pin eat ad 4 B ¥ - M eee 7h ae areete sk tyret Ld 

f : , renee ‘ CP as Ay Pas td he - i Pi heh eat or 

Ret ; chee teh 
of 

yah 


Pie re wae 
i 
Ukiah 
i 


ML $4 < by 
ESRD SS ere reise ke) LMedet HK ghee fe ey 
Uy bt hoe eh wee rad v | omg ts oy is s 
ne ea ae Rane Hy Seite rie hr Peas age te Seeera mene 
ga 1 Ras pei Eee Pa ietet ht pe kh tba or PS 4 
REC NS tal 2 faci . eee Se PPE eta haat Bet Lads Lah reife 
pet Melt fata 9 Me peti Faas | ek Lek bih Lab pee ABO leh 4 eae ats a 
OM REMC TREY Feed tal Tpegesepearesesst ere: 
ARC Re CR AK ott CP a ita et 
GB fs f ‘J a 4 a ct Wet a 
Rt ae SE Pra Cais te Leta 
Ot ete) cr) aes Vata ces 20s oe My 
En Coles Stas ih ae At Cia ie 
Tet Pras ipa OT be 
Paros esau hy A Lees t rare 
fay ee ae el ied | _ 
A F he 
‘ 


it 
Bre 
Pe 


ie 
pee 
if tae- es 

Laas Js 


ea 
PT tr hehe by tt ti 
bi, 9 





DUDLE’ “NOX LIBRARY 
NAVA) 3 TGRADUATE SCHOO! 
MONTEREY GA 93943-5101 








Form Approved 
REPORT DOCUMENTATION PAGE ae ee 


Public reporting burden for this collection of information ts estimated to average | hour per response, inctuding the time for reviewing instructions. searching @risting data sources, 
gathering and maintaining the data needed, and completing and reviewing the cottection of information Send comments regarding this burden estimate or any other aspect of this 


collection of information, including suggestions for reducing this Durden to Washington Headquarters Services, Directorate for information Operations and Reports, 1215 Jefferson 
Davis Highway, Suite 1204, Arlington, VA 22202-4302 and to the Office of Management and Budget. Paperwork Reduction Project (0704-0188). Washington, OC 20503 
1. AGENCY USE ONLY (Leave blank) 2. REPORT DATE 3. REPORT TYPE AND DATES COVERED 

Mange 993 Master's Thesis 


4. TITLE AND SUBTITLE 5. FUNDING NUMBERS 
Resolution in Radar Mapping 


6. AUTHOR(S) 


Michael D. Anderson 





7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 
Naval Postgraduate School 
Monterey, CA 93943-5000 


8. PERFORMING ORGANIZATION 
REPORT NUMBER 















9. SPONSORING /MONITORING AGENCY NAME(S) AND ADDRESS(ES) 10. SPONSORING / MONITORING 


AGENCY REPORT NUMBER 


Pe eeeM ery NO’ The viewS expressed in this thesis are those of the 
author and do not reflect the official policy or position of the 
Department of Defense or the US Government. 


12a. DISTRIBUTION / AVAILABILITY STATEMENT 12b. DISTRIBUTION CODE 






















Approved for public release; distribution is 
unlimited 








13. ABSTRACT (Maximum 200words) Signal processing has led to great performance 
gains in radar mapping. The most critical feature of these systems is 

cell size, which determines resolution. Cell size is defined by range 

resolution and azimuth resolution. 

Range resolution is improved through pulse compression. Phase of 
frequency modulation of a waveform yields increased bandwidth and short 
er effective pulse width without reducing total signal energy. Several 
families of codes are investigated emphasizing matched filter output 
and doppler tolerance. 

Azimuth resolution is improved through beam sharpening. Several 
beam sharpening techniques are illustrated with radar images provided 
by Hughes Aircraft. Range bin output plots demonstrate the effective- 
ness of these techniques. 

With these techniques, "near-SAR" quality output can be obtained 
from real beam mapping radars allowing the real-time and all aspect 
capabilities of real beam systems to be more fully employed in 


tactical missions. 
15. NUMBER OF PAGES 
90 


14. SUBJECT TERMS 
16. PRICE CODE 






Beam Sharpening; 











Resolution; Pulse Compression; 
Radar Mapping 









18. SECURITY CLASSIFICATION [| 19. SECURITY CLASSIFICATION {| 20. LIMITATION OF ABSTRACT 


OF THIS PAGE OF ABSTRACT 
UNCLAS UNCLAS UL 


NSN 7540-01-280-5500 Standard Form 298 (Rev 2-89) 
2 eapreslriBed go" 4 NSTSta~239 18), 





17. SECURITY CLASSIFICATION 
OF REPORT 








UNCLAS 





Approved for public release; distribution is unlimited 
RESOLUTION IN RADAR MAPPING 
by 
Michael D. Anderson 
Lieutenant, United States Navy 
B.S.E.E., United States Naval Academy, 1987 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN ELECTRICAL ENGINEERING 


from the 


NAVAL POSTGRADUATE SCHOOL 
March 1993 


om wD? 


ae 


ABSTRACT 

Signal processing has led to great performance gains in radar mapping. The 
most critical feature of these systems is cell size, which determines resolution. Cell 
size 1s defined by range resolution and azimuth resolution. 

Range resolution is improved through pulse compression. Phase or frequency 
modulation of a waveform yields increased bandwidth and shorter effective pulse 
width without reducing total signal energy. Several families of codes are investigated 
emphasizing matched filter output and doppler tolerance. 

Azimuth resolution is improved through beam sharpening. Several beam 
sharpening techniques are illustrated with radar images provided by Hughes Aircraft. 
Range bin output plots demonstrate the effectiveness of these techniques. 

With these techniques, "near-SAR" quality output can be obtained from real 
beam mapping radars allowing the real-time and all aspect capabilities of real beam 


systems to be more fully employed in tactical missions. 
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I. INTRODUCTION 

The evolution of radar has reached a new stage. Like many other areas of technical 
development, software is now driving innovation and providing new levels of performance. Digital 
signal processing has permeated every aspect of radar, from the shape of the transmitted pulse to 
the type of return display. One application which has greatly benefitted from signal processing 
advances is radar mapping. 

A. RADAR MAPPING 

One of the primary advantages of radar mapping 1s all-weather capability. Clouds, fog and haze 
generally do not limit system capability. In addition, radar mapping can be done day or night. 
Finally, both the angle and direction of illumination are controllable. For these reasons, radar 
mapping is invaluable in environmental mapping as well as tactical surveillance. Airborne radars are 
able to fill a number of mission roles through mapping. 

Navigation is one of the roles for radar mapping. Numerous applications exist to fill this role. 
Terrain following involves vertical scanning with horizon information fed directly to flight control to 
allow extreme low altitude approach. Terrain avoidance is similar, but incorporates horizontal 
sweeps to avoid obstacles. TERCOM, or terrain contour mapping, is used to guide a platform over 
a timed preprogrammed ground trajectory. This mode assumes a known contour map. Radar map- 
matching involves guiding a platform over a pre-designated course over previously mapped terrain, 
as the Tomahawk TLAM cruise missile does. Either real-beam or SAR mapping can be used for 
this function. 

In air-to-ground applications, weapon delivery can be entirely dependent on radar mapping. All 
weather attack capability is demonstrated in the scenario of Figure 1. [1] In this case, a SAR 
processing system is demonstrated. A first SAR map gives the system operator an overview of the 


attack area. By designating a point with a cursor, a second finer-detail image is mapped of the 


1 


immediate target area. By marking the precise target with a cursor, the operator initializes a 
sequence in which the system provides steering guidance to the target and controls weapon release. 

In the air-to-ground battle, radar mapping also proves invaluable in reconnaissance and 
surveillance. The ability to distinguish and localize targets from stand-off distance provides a great 
advantage to friendly ground forces. Long-range targeting data and advance warning of large troop 
movements can provide the margin of victory in battle. Finally, radar mapping can also be quite 
useful in the area of battle damage assessment. 

For several of these mission roles, SAR mapping has been the method of choice. With recent 
developments in real-beam mapping, however, real beam is expanding its role in warfare. The value 
of real beam mapping was demonstrated in Operation Desert Storm. A vital mission of the U.S. Air 
Force and other allied air assets was the location and destruction of both fixed and mobile ground- 
to-ground missile or Scud launchers. This became a high-priority cat and mouse game as efforts to 
protect Israel and Saudi Arabia from missile attack required searching great expanses of desert for 
relatively small targets. In effect, this mission resembled the proverbial search for a needle in a 
haystack. While Synthetic Aperature Radar (SAR) could provide the necessary resolution to detect 
these launchers, the time required for processing of the data would slow the searches considerably. 
With its all-aspect real time capability, real beam systems could search a much greater area in less 


time. This is critical when dealing with mobile targets and quick response tactics. 
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As envisioned, a real beam system would perform rapid sweeps of the desert and localize 
potential targets and areas of interest. A SAR system would then be employed to provide greater 
detail and targeting information. In this way, limited assets could be best employed to meet the 


threat as quickly as possible. 


B. CELL SIZE AND RESOLUTION 

As alluded to above, the effectiveness of radar mapping 1s, to a large extent, dependent on 
system cell size at operational range. Cell or "pixel" size is defined as a rectangle whose sides are 
the system range and azimuth resolution distances. Driving the choice of cell size are considerations 
such as object size, processing requirements, interpretation requirements and, of course, cost. A 
demonstration of the effects of cell size is presented in Figures 2 and 3. {1] In this figure, all 
portions of the shape were assumed to reflect equally with full cells lighter in the corresponding map 
and partially full cells darker according to the strength of the return. In actuality, several mappings 
from different aspects with different frequencies and polarizations would have to be combined to 
generate this quality return. Obviously, with smaller cel! size, detection and localization of small 
objects is improved as well as discrimination of detail in larger objects or areas. Table 1 [1] provides 
typical cell size requirements for square mapping of common map features. 

In general, it is easier to achieve fine range resolution than fine azimuth resolution. Resolution 
in range is approximately 492 feet per microsecond of pulse width. In order to improve range 
resolution beyond this figure, all that needs to be done is to reduce the pulse width of the system. 
As will be discussed in Chapter II, when there are limits to the physical reduction of pulse width, 
pulse compression can be used to achieve additional improvement in resolution. The matched filter 
and ambiguity diagram, as presented and implemented in Appendices A and B, will be used to 
illustrate several techniques. However, the more significant limitation for real beam systems has 


always been the azimuth resolution. 


CORRESPONDING MAP 
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CELL SIZE’ 1/5 MAJOR DIMENSION 


Figure 2. Simulated resolution of cell size equal to 1/5 object size. 
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Simulated resolution of cell size equal to 1/20 object size. 


Figure 3. 


TABLE 1 
CELL SIZE REQUIREMENTS 


Features to be Resolved Cell Size 


Coastlines, large cities, outlines of mountains 
Major highways, variations 1n fields 60-100 ft 


"Road map” details: city streets, large buildings, 30-50 ft 
small airfields 


Vehicles, houses, small buildings 5-10 ft 


Azimuth resolution of a standard system is roughly the 3 dB beamwidth of the antenna times 













the range. At five nautical miles, a 16 foot array operating in X band has an azimuth resolution of 
approximately 200 feet. To recognize the shape of even large objects at long distance, the military 
has had to resort to SAR. Since the beamwidth varies proportionally with wavelength and inversely 
with antenna length, to improve resolution, it is necessary to operate at shorter wavelengths or with 
larger antennas. There are limits on choice of wavelength and antenna size. Limitations on 
machining tolerance and phase alignment make the construction of large antennas a difficult task. 
Airframe requirements complicate this task even further. Improvements in signal processing, 
however, have advanced to the point that "near SAR" quality images can be obtained from real beam 
systems. This is especially true for point reflectors. Chapter III will present a discussion of azimuth 
resolution and several beam sharpening techniques. These techniques will include actual beam 
sharpening algorithms as well as image enhancement techniques which are applied after initial beam 
sharpening. The comparative performance of these techniques will be demonstrated by the quality 
of through a set of images obtained through the application of algorithms developed by the 
Advanced Systems Division of Hughes Aircraft Company. 

As these techniques are further developed and refined, real beam radar mapping will 


undoubtedly be applied to a great many situations which require SAR systems today. With faster 


scan time, all aspect capability and good resolution, real beam systems can provide a great advantage 


in tactical operations. 


II, RANGE RESOLUTION 


A. WAVEFORM SELECTION 

Long operational range and fine range resolution are two of the primary design considerations 
in a mapping radar. These apparently contradictory goals are achieved by proper selection of radar 
waveform. The radar waveform is defined by its carrier frequency, pulse repetition frequency 
(PRF), pulse width and intrapulse modulation. The tradeoffs involved in the selection of waveform 
parameters are described in this section. 

The selection of carrier frequency is an important first step in waveform design. The 
propagation characteristics of radar change rapidly in the microwave region. Nathanson [2] provides 
the following guidance. Attenuation due to rain (in dB) is proportional to f?® and backscatter from 
rain and small particles varies as {* dB over most of this region. Ionospheric effects also vary with 
frequency and are of concern at frequencies below 3 GHz, while polar operations are affected by 
aurora backscatter below 2 GHz. Table 2 [1] defines the radar frequency bands and common 
usages. The carrier frequency also has implications for azimuth resolution, though this will be 
discussed at greater length in the next chapter. 

Once the operating frequency has been chosen, the next issue to be considered is the selection 
of pulse repetition frequency (PRF). The most widely used method of range measurement is pulse 
delay ranging. It is simple, avoids the problem of transmitter leakage noise, provides easy range 
measurement, and is extremely accurate. With this form of operation, however, great care must be 
taken when selecting the PRF. The choice of PRF is crucial because it determines to what extent 
the ranges and doppler frequencies observed by the radar will be ambiguous. 

The returns from objects separated by the system unambiguous range are received 


simultaneously. The echoes from a target, therefore, must compete with ground return not only 


from the target's own range but from every range that is separated from it by a whole multiple of 
the unambiguous range. Any overlapping of the return profiles can result in target echoes and 
ground clutter passing through the same doppler filter, even though the true doppler frequencies of 
the target and clutter are quite different. This becomes a significant problem if the received 
frequency is less than the width of the true doppler profile, as it must often be to reduce range 
ambiguities. Reducing the PRF causes the mainlobe clutter to occupy a large portion of the receiver 
passband. As the percentage of the passband occupied by the mainlobe clutter increases, it becomes 
difficult to reject clutter on the basis of doppler frequency without rejecting a large percentage of the 
target echoes as well. The lower the PRF, the more severe this effect becomes. With ground 
clutter being widely dispersed, range and doppler ambiguities greatly compound the problem of 
isolating target echoes. 

Once the decision regarding PRF has been resolved, the designer's attention turns to the 
selection of pulse width. Ideally, to achieve both long detection range and fine range resolution, a 
system should transmit extremely narrow pulses of exceptionally high peak power. Since there are 
limitations on the amount of peak power which can be transmitted, to obtain long range at PRF's 
low enough for pulse delay ranging, fairly wide pulses must be transmitted. Pulse length determines 
the range resolution. To improve range resolution, a designer must shorten the pulse width. As 
pulse length is decreased, however, so is the total energy contained in the pulse. Eventually, a point 
is reached where further energy decrease in the signal is unacceptable due to detection 
requirements. With transmitter peak power and receiver sensitivity improvements coming slowly, 
and only at great cost, a clear limit seemed to exist in the detection-range resolution tradeoff. In 
search of a way around these limits, it was found that the limits on range resolution can be 
circumvented by coding successive increments of the transmitted pulse with phase or frequency 


modulation. This led to the development of linear frequency-modulation pulse compression. 
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TABLE 2 
RADAR FREQUENCY APPLICATIONS 


Bands Frequency Usage Summaries 
Range 
3-30 MHZ Over the horizon (OTH) Radar. Very-long-range with lower 
resolution and accuracy. 
30-300 MHZ Long-range line-of-sight (LOS) surveillance. Low-medium 
UHF 300-1000 MHZ resolution and accuracy. No weather effects. 


HF 
VHF 
Cc 1-2 GHZ Long-range surveillance/enroute traffic control. Medium 
resolution, slight weather effects. (200nmi) 
S 
Cc 
Xx 
K, 
K 















1-2 GHZ Long-range tracking w/high accuracy. Airborne weather 
detection. Increased weather effects in rain. 
8-12 GHZ Short-range track. Missile Guidance. Mapping/weather 


12-18 GHZ High-resolution/synthetic aperture mapping. Satellite altimetry 
Short-range or good weather required. 


18-27 GHZ Limited use due to water vapor attenuation. 


Very-high-resolution mapping. Airport surveillance. 
40-300 GHZ Smart munitions, remote sensing. Weather limited. 


2-4 GHZ Short-range surveillance. Terminal air traffic control. Moderate 
weather effects. (60nm1) 
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If a matched filter is used, the sensitivity of a radar receiver depends only on the total energy 
contained in the signal. The time distribution of the energy is unimportant. Within the constraints 
set by frequency and PRF selection, resolution then becomes the determining factor in the form of 
the transmitted energy. The single measure which best expresses the range resolution properties of 
a signal is the effective bandwidth. Long duration and high bandwidth are not incompatible if a 
signal contains rapid or irregular structure changes such as in pulse compression. Three types of 
pulse compression will be discussed in this chapter, linear frequency modulation, pseudo-random 


binary coding and pseudo-random polyphase coding. 


B. LINEAR FREQUENCY MODULATION 

Linear frequency modulation (LFM), also known as "chirp," 1s the result of increasing or 
decreasing the radio frequency of a transmitted signal at a constant rate throughout each pulse. This 
form of pulse compression enables very large compression ratios, and is simple to implement. There 
is a slight ambiguity between range and doppler using this method. A positive doppler shift will 
cause the signal to leave the filter sooner. The radar system can not distinguish between a target at 
a slightly closer range and one with a positive doppler shift at a slightly farther range. 

To resolve a LFM signal, it is passed through a linear filter with a frequency change equal to the 
inverse of the pulse width. This must occur over a time equal to the compressed pulse width. The 
pulse compression ratio (PCR) achieved is equal to the ratio of the transmitted pulse width to the 
compressed pulse width. This is also equal to the pulse width times the frequency change which is 


known as the "time-bandwidth product." 


C. PSEUDO-RANDOM BINARY CODES 
In many applications, digital implementation of pulse compression is preferred over analog. 


There are many advantages to digital implementation, including: 
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- freedom from ringing due to impedance mismatches/analog filters 
- reproducible response 

- high peak-to-sidelobe ratios without weighting 

- change bandwidth by changing clock frequency 


- change the waveform and PCR by controlling the number of time samples and by digitally 
switching the matched filter 


- compatibility with MTI and other signal processing. 

Despite the advantages of digital implementation, significant limitations still exist due to range 
sidelobes. Range sidelobes exist due to the sin(x)/x shape of the energy spectrum of the 
uncompressed pulse. As the uncompressed pulse passes through the filter, the higher frequency 
spectral sidelobes travel faster than the main lobe frequencies, and the lower frequency sidelobes 
travel slower. The spectral shape of the energy then is reflected in the shape of the compressed 
pulse amplitude vs time plot. This has an effect similar to antenna sidelobes. Range sidelobes imply 
wasted energy and can result in false targets. Amplitude weighting, such as the use of Hamming 
window, can be used to minimize the sidelobes. This weighting, however, often has the effect of 
broadening the compressed pulse width or “main lobe" of the energy form. 

The binary coding is actually a phase reversal of 0 or 180 degrees impressed upon the carrier 
signal. These phase sequences are usually represented by a series of +1 and -1. The resulting 
signal is wideband with bandwidth equivalent to twice the code clock rate. The range resolution of 
these codes is determined by the transmitted chip width. A chip is defined as one coded segment of 
the transmitted pulse. 

The following are several of the more popular binary phase coding sequences. Most are easily 
implemented and display good autocorrelation characteristics. This class of waveform, however, 


does suffer from significant doppler phase shift limitations. 


1. Barker Codes 

These are a select number of pseudo-random codes which come close to the desirable feature 
of zero range sidelobes. With these sequences, the sidelobes of the compressed pulse are limited to 
the amplitude of the uncompressed individual segment, while the main lobe amplitude is multiplied 
by the number of bits in the code. Matched filter outputs for these sequences have normalized 
range sidelobes of only a/li2 relative to the mainlobe. The primary drawback of these codes is that 
the longest known code is only 13 bits long. It has been proven that no other examples of these 
codes exist for odd lengths, and no even codes for lengths less than 6084. While it has not been 
proven, it is probable that no others exist at all. The known Barker codes, along with their known 
autocorrelation peak sidelobes and integrated sidelobes are given in Table 3. [2] The limitation to 
length 13 is serious, since it does not allow complete decoupling between average power and 
resolution. A second drawback is that the matched filter output of these waveforms degrades rapidly 
with doppler shift. Large sidelobes result in the autocorrelation and appear as spikes other than at 
the origin in the ambiguity diagram. This effect is demonstrated in Figure 4 which shows the 
matched filter output without doppler shift and the ambiguity diagram for a Barker 13 sequence. 


TABLE 3 
THE KNOWN BARKER CODES 


Code Elements PSL, dB | ISL, dB 
as PS a 


pos teers | sso | so 


ee ee -20.8 -10.8 


a 
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2. Compound Barker Codes 
One attempt to circumvent the length limitation of the Barker codes is known as compound 


Barker codes. This coding process can best be described through the use of the example below: 


Barker length 4: + + - + — sidelobe level -12.0dB 
Barker length 3: => (++-) (++-) (++-) (++-) sidelobe level - 9.5dB 
Compound Barker length 12; ++- ++- +--+ ++- sidelobe level - 9.5dB 


Each bit of one sequence represents either a true or inverted replication of the second sequence. 
The resulting length and pulse compression ratio is equal to the product of the two original code 
lengths. While the range sidelobes will no longer be equal in amplitude, the main lobe to sidelobe 
amplitude ratio will be equal to the lowest original code used. Figure 5 provides the matched filter 
output and ambiguity diagram for Complementary Barker length 25 (5 in 5). This figure 


demonstrates that these longer codes are very sensitive to doppler phase shift. 
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Mreure 5. 
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3. Minimum Peak Sidelobe (MPS) Codes 

Another effort has been the compilation of the minimum sidelobe codes. These are simply 
the binary sequences which attain the lowest PSL for given lengths. Table 4 [2] gives the MPS codes 
for lengths 14 through 48. These codes were discovered by an exhaustive search of all possibilities. 
Figure 6 gives the matched filter outputs and ambiguity diagrams for the MPS code length 5 (a 
Barker) and a nonideal code of the same length for comparison. Note the higher peak sidelobes in 
the non-optimal code. 

4. Complementary Codes 
There are four known kernels of complementary code. A kernel is a basic length of code 
which cannot be shortened by inverting standard generation procedures. The known kernels are 
given in Table 5 [3] below. Corresponding sidelobes produced by the A and B sequences have 
opposite phases. By alternately modulating successive pulses and switching the matched filter, 
sidelobes are canceled upon integration of successive pulses. 

Longer sequences can be built up by “chaining” to form what is known as composite 
complementary sequences. Chained codes have sidelobes greater than a single code, but these 
cancel upon integration of an even numbers of pulses. The topic of complementary codes and their 
ambiguity diagrams were explored in Akita [3]. Figures 7 and 8, from this source, present an 
example of matched filter output and the corresponding ambiguity diagram for a complementary 


sequence of length 10. 
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TABLE 4 
MINIMUM PEAK SIDELOBE (MPS) CODES 








Sample code 
Length Number PSL ISL, dB (octal) 
a ] 1 -—9.12 047 
8 16 2 - 6.02 227 
9 20 2 — §.28 327 
10 10 2 — 5.85 0547 
1] 1 1 — 10.83 1107 
12 a2 2 — 8.57 4657 
13 M 1 -11.49 12637 
14 18 2 -7.12 12203 
15 26 2 - 6.89 14053 
16 20 2 ~ 6.60 064167 
1 8 Z — 6.55 073513 
18 4 2 -8.12 310365 
19 2 2 — 6.88 1335617 
20 6 Z = 121 1214033 
21 6 2 - 8.12 5535603 
22 756 3 — 7.93 03466537 
23 1021 3 — 7.50 16176511 
24 1716 3 — 9.03 31127743 
25 2 2 —8.51 111240347 
26 484 3 — 8.76 216005331 
27 774 3 — 9.93 226735607 
28 4 pi - 8.94 1074210455 
29 561 3 = 6.31) 2622500347 
30 12 3 — 8.82 4305222017 
31 502 3 — 8.56 05222306017 
32 844 3 — 8.52 00171325314 
33 278 3 — 9.30 31452454177 
34 102 3 -— 9.49 146377415125 
35 222 3 - 8.79 000745525463 
36 322 3 — 8.38 146122404076 
a7 110 3 — 8.44 0256411667636 
38 34 3 =9 19 0007415125146 
39 60 3 — 8.06 1146502767474 
40 114 3 - 8.70 02104367035132 
41 30 3 — 8.75 03435224401544 
42 8 3 - 9.41 04210756072264 
43 24 3 — 8.29 000266253147034 
he 30 3 — 7.98 017731662625327 
45 8 3 =o 052741461555766 
46 2 3 — 8.12 0074031736662526 
47 2 3 — 8.53 0151517641214610 
48 8 3 sate | 0526554171447763 








a) Ambiguity diagram for the MPS length 5 code (+++-+). 
b) Ambiguity diagram for a non-optimal code of length 


Figure 6. 


5 (+4+--+). 
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TABLE 5 
THE KNOWN COMPLEMENTARY KERNELS 


LENGTH | SEQUENCE 


= {-1,+1,+1,-1,+1,-1,+1,+1,+1,-1} 
= {-1,+1,+1,+1,4+1,+1,+1,-1,-1,+1} 


Aare 1-1 leaigs 1+ 1 ta ey 
B= 4+1541,4+1,+1-1,+ 1, 41-l el} 


=a F 1 1)- eae ieee 
+ 1,-1,-1,-1,-1, +1, + 1,-1,-1,-1, +1,-1,+1} 

B = {-1,+1,-1,-1,+1,+1,-1,+1,+1,+1,+1,-1-1, 
-1,-1,-1,-1,-1,+ 1,+ 1,-1,-1,-1,+1,-1,+1} 





The complementary codes, like all binary phase codes, are limited by doppler sensitivity. Any 
doppler shifts encountered must be small, or pulse widths must be short, since it is the total phase 
shift across the pulse which matters. Excessive phase shift across the waveform will result in serious 


degradation of the amplitude vs. time plot. 
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combined output 





Figure 7. Matched filter output of complementary length 10 kernels; 
sequence A, sequence B and summation. 
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Figure 8. Ambiguity diagram of complementary length 10 kernel. 
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C. POLYPHASE CODES 
Polyphase codes are another class of pseudo-random sequences. With these codes, the bits can 
consist of any number of different harmonically related phases, rather than only 0 or 180 degrees. 
Several such codes are available, with some of the more common ones described below. 
1. Frank Codes 

The Frank codes are a series of inphase and quadrature (IQ) samples taken at the Nyquist 
rate of a “step-chirp" waveform. The formation of the code is demonstrated by Figure 9. [4] To 
generate this code, P phase increments are defined by dividing 360/P. P groups of P segments are 
then generated according to the following procedure: 

The first phase of each group is 0. The phases of the remaining segments in each group 


increase in increments of 


A® =(G-1)*(P-1)*8 (2.1) 
where 
G is the group number and 


@ is the basic increment. 


The resulting code is given by 


@ (i,j) = (2x /N)(i-1)(-1) (2.2) 
For example: 
P=4; 
Group 1: segment 1=0, segment 2=0, segment 3=0, segment 4=0 
Group 2: segment 1=0, segment 2=(360/4)(1)(1) =90, segment 
3=(360/4)(1)(2) = 180, etc. 


The resulting code is shown at the bottom of Figure 9. 
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Figure 9. Derivation of Frank code from a step-chirp waveform. 


ZS 


For a given number of segments, a Frank code provides the same pulse compression ratio as 
a binary phase code and the same peak to sidelobe ratio as a Barker code. Yet, by using more 
phases (increasing P), the codes can be made any length. As P is increased, however, the size of the 
fundamental phase increment decreases, making performance more sensitive to externally introduced 
phase shifts and imposing more severe restrictions on uncompressed pulse width and maximum 
doppler shift. This is demonstrated in Figures 10 and 11 showing degradation of the ambiguity 
diagrams for increasing length Frank codes. Another limitation of the Frank codes is that the 
largest phase change between elements occurs at the center of the sequence. Because of this, 
receiver bandwidth limitations have especially significant effects on the center peak to sidelobe ratio. 
2. P1 Code 

To alleviate the problems of bandwidth limitation present in the Frank codes, the code groups 
representing the } frequencies were rearranged in time order of transmission. This placed the lowest 
amount of phase increment in the center to preserve the main lobe. This led to the creation of the 
P1 code given by 

@ (ij) = -(e /N)[N-(2-1)]1G-I)N+ FD) (2.3) 
When N is odd, the P1 code is simply the Frank code rearranged to have conjugate symmetry 
around the DC term. The autocorrelation functions of the two codes are nearly identical, but the P1 
has no bandwidth limitations. 
3. P2 Code 
The P2 code is given by 
@ (ij) = {(% /2)[(N-1)/N]}-(x /N)(i-1)} {N+1-2)} (2.4) 

It has an autocorrelation function nearly identical in magnitude to both the Frank and P1 codes, and 


has no bandwidth limitations. The difference is that P2 is real as opposed to complex 





(length 9) and p=4 (length 16). 


Figure 10. Frank codes of P=3 
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Figure ll. Frank codes of P=5 (length 25) and P=7 (length 49). 
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because it is symmetrical. Figure 12 [4] demonstrates the improved bandwidth tolerance of both P1 
and P2 over the Frank codes. In this example, bandlimited results show a 6 db improvement in peak 
sidelobe level for P1 over Frank code. Both P1 and P2 still remain fairly doppler intolerant, 
however. This is a characteristic of the analog step chirp waveform from which both, along with the 
Frank code, are derived. 
4, P3 Code 
Since the linear chirp waveform is much more tolerant of doppler than the step chirp, the P3 
code was developed from it. The code is written as 
@(i) = x(i-1)”/p as (i) varies from 1 to p (2.5) 
In this equation p represents the pulse compression ratio. The resulting waveform is just as doppler 
tolerant as the analog linear chirp, but suffers from the same bandwidth limitations as the Frank 
codes. Figure 13 demonstrates the improved doppler tolerances of the P3 code. While the unshifted 
matched filter output shows the Frank to have a superior peak-to-sidelobe ratio by 4 db, the doppler 
shifted signals clearly show the superiority of P3. Of particular note is the elimination of grating 
lobes and preservation of the image lobe. 
5. P4 Code 
A waveform both doppler tolerant and not bandwidth limited resulted from rearranging the P3 
code. The resulting code is known as P4 and is written as @(1) = [x (i-1)7/p ]-t (i-1). Figure 14 [4] 


demonstrates the bandwidth limitation improvement of P4 over P3. 
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Figure 12. Effect of precompression bandlimiting of Frank and Pl codes 
of length 100. 
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Figure 13. a) and b) are Frank and P3 code matched filter outputs for 
length 100. c) and d) show matched filter outputs for the 


doppler phase shifted case. 
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Figure 14. Bandlimited matched filter outputs for P3 and P4 codes length 100. 
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II. AZIMUTH RESOLUTION 

With pulse compression providing a means of achieving good detection range and fine range 
resolution, the largest remaining constraint in a mapping radar ts azimuth resolution. As mentioned 
in the introduction, in an unprocessed real beam system, azimuth resolution is determined by the 
antenna beamwidth. While beamwidth control is adequate for target detection, prohibitively large 
antennas result for the narrow beams required to determine detailed characteristics such as the 
number, size or separation of targets. The need for information of this sort led to the development 
of synthetic aperature radar. 

Unfortunately SAR systems also have limitations. These systems are unable to map at nose 
aspect, are expensive, and are hardware dependent, Most importantly, they are comparatively slow 
in mapping due to the extensive amount of signal processing required. These limitations have led to 
renewed efforts to improve the angular resolution of real beam systems. These efforts are referred 


to as real beam sharpening. 


A. BEAM SHARPENING TECHNIQUES 

Beam sharpening reflects two areas of effort: actual beam sharpening algorithms, and image 
enhancement techniques which are applied after initial beam sharpening. While some 
implementations will combine some of both, they will be presented individually in this paper. 

These techniques enhance angular designation accuracy in all aspect scenarios. Further, real- 
time post-processing is possible to accomplish this without extending the radar frame time. Even 
while implementing these additional processors, real beam systems can map a 60 degree swath in 
just 1-2 seconds. 

Techniques already in use on the F-15E provide a 2-1 enhancement in designation accuracy. 


Test results of new techniques being developed at Hughes Aircraft promise beam sharpening of 6:1 
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with 6 dB SNR and up to 25:1 for a 20 dB SNR return. These estimates are based on isolated 
corner reflector return. The following pages will provide a brief introduction into a number of these 
techniques. 

1. Combination of Sum and Difference Channels 

Useful for resolving point targets, this method was developed at Hughes Aircraft and is 

currently implemented in the F-15E. Figure 15 [5] demonstrates the concept. The boresight null in 
the difference (D) channel pattern forms a sharp peak when combined with the sum (S) channel 
pattern. A weighting factor (K) controls the beam sharpening ratios. A value of (K) equal to .75 is 
common. The output (S,) is given by 


s. = |S] - KID] (3.1) 


This technique can also be applied using sum and synthetic difference channels. If difference 
channel information is not available, synthetic difference channel data may be generated through the 
use of anti-symmetric weights. The difference channel data (D) is generated as shown below. 

Be 2X (1-k) w, Si+k (3.2) 
where w, is a Suitable set of anti-symmetric weights, (j) is the radar pulse index, and oer is the sum 
channel data for each pulse and set of weights. The optimum weights (w,) to maximize the gain 
slope of the synthetic difference pattern vary as 

w,« dG@)/dé (3.3) 
where: 

G(@®) = antenna azimuth gain 

@ = k(scan rate)/PRF 


The sum and synthetic difference channel data are then combined as for the previous case. 


5 
= 
V) 
2 
o 
ce 
8 
= 
& 
V) 
2 


a. Difierence and Sum Channel Magnitude 





Figure 15. Combination of sum and difference channels resulting in 
attenuation of response away from boresight. 


5D 


2. Monopulse Beam Sharpening 

Also used for point targets, this technique uses sum and difference channels to form an 
azimuth discriminant on the I/Q data for each return in each range bin. The discriminant (a) for 
each pulse (j) is given by 

a= kRe((S; .D,)/(S, S)) (3.4) 

where k is a scaling factor given by 1/(discriminant slope)(angular bin spacing), Sj is the sum 
channel data for each pulse and D; is the difference channel data for each pulse. Re is the real 
operator and * denotes complex conjugation. This discriminant value gives the distance of each 
scatterer from the center of the range bin. As the scatterer enters the beam, each subsequent pulse 
will show the scatterer moving closer to the center as the beam scans across it. Eventually, the 
discriminant will reach zero, cause a sign change and move out of the beam. Once the discriminant 
is formed for a pulse and range bin, the value of the discriminant assigns the return to one of the 
azimuth bins which each range bin 1s divided into. The magnitude of the sum channel data is what 
is actually stored in the bin. With each pulse, data is continually added to the bins. Upon 
completion of the scan, the accumulated sum of the contents of each bin is displayed thereby 
localizing scatterers. 

Synthetic monopulse beam sharpening ts also possible. It involves the sum channel and 
synthetic difference channel. The synthetic difference channel is formed as previously described, and 
is then applied as in regular monopulse beam sharpening. 

3. Inverse Filtering 
In a large class of imaging systems, spatial degradation can be modeled by a linear-shift- 
invariant impulse response and additive noise. In these cases, restoration can be accomplished 
through linear filtering. This is illustrated as part of an algorithm in Figure 16. [6] The image 


degradation is modeled by a system with a particular transfer function Hp(x,y). Noise N(xy) is 





(Ax) By (Ain) Fy 
NOILvaVYD30 


JDVWI 
qVivas 





(txt 





(Any ty 





Y3ligs 
NOILVYOLSSY 


a 


dt 2S Ge (NS @ (A ak aR) aA er) Oe 


JO 


(AX) + gpep(g —-dx—x)"y(g-x)ia || = (Axe 
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added to the signal resulting in the observed image field Fo(xy). After passing through a 
restoration filter, an approximation to the original image field F is obtained. 

Inverse filtering reflects the earliest attempts towards image restoration. In this process, the 
transfer function of the degrading system is inverted to yield a restored image. One application of 
this technique is the modeling of low-resolution data as the convolution of high resolution data with 
a known spatial antenna pattern. With a known transfer function for the antenna and a good 
estimate of the mean square of the observed data, the mean square of the high-resolution signal can 
be reconstructed through appropriate deconvolution techniques. This assumes that the original and 
observed data are random processes. If source noise is present, however, an additive reconstruction 
error will result. This error can be quite large at spatial frequencies for which the spatial image 
degradation transfer function is small (typically, at high spatial frequencies). This impairs the quality 
of high detail regions of an image. Proposals have been made to reduce susceptibility to noise. One 
method involves using a restoration filter with a transfer function which is multiplied by a unit step 
function over the areas where the image spectrum is expected to exceed the noise spectrum, and ts 
zero everywhere else. This strikes a compromise between noise suppression and loss of high- 
frequency image detail. Other filter implementations, such as the Wiener filter and its variations, 
provide better performance in the presence of noise. 

4. Maximum Entropy 
Maximum entropy is a nonlinear technique which can be used to resolve closely spaced signals. 
It is an extension of inverse filtering which incorporates linear prediction to better resolve an image 
scene. The technique is applied to magnitude data in each range bin. As scatterers enter the beam, 
they cause predictable transients due to the roll-off of the antenna beam shape. A Fourier 
transform of the data is then taken and scaled to remove the effect of the antenna pattern. This 
results in a set of autocorrelation coefficients. A set of predictive linear equations ts applied and 


solved to further resolve the data which has been corrected for actual antenna length and generate a 
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second set of coefficients. A second Fourier transform is then performed, this time on the 
coefficients. The resulting transform is related to image intensity with the poles of the coefficients 
denoting the scatterers. 

5. Hughes Advanced Discrimination Technique (HADT) 

Further work has been done in this field by personnel at Hughes Aircraft Advanced Systems 
Group. While much of their work is proprietary in nature and cannot be discussed in detail, the 
following figures show the improvement in resolution which can be achieved with today’s signal 
processing. 

Figure 17 [5] provides a comparison of the various beam sharpening techniques for both point 
and extended targets. Note that for point targets, the maximum entropy, monopulse beam 
sharpening and HADT methods provide excellent results. For extended targets, the maximum 
entropy and HADT work well with the exception of some peaking in the return just before roll-off. 
A Hughes modified inverse filter has also attained good results on extended targets. (Figure 18 [5] 
shows the real return of an air-to-ground radar with a 2.5 degree beamwidth mapping two corner 
reflectors spaced three degrees apart.) Even though the separation is greater than the 3 dB beam 
width, the objects are indistinguishable in the raw return. In the HADT processed response, 
however, the targets are readily distinguished. Figure 19 [5] provides a comparison of SAR results, 
unprocessed real beam return and sharpened real beam ground mapping (RBGM) return. While 
targets in the unprocessed real beam image appear as smears due to the wide beam width sweeping 


across them, the HADT sharpened image reveals the SAR processed image in clarity. 
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Figure 17. Relative effectiveness of commom beam sharpening techniques on 
point reflectors and extended objects. 
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and 





SAR Map (Resolution 42 ft) 


RBGM After Beam Sharpening 








B. IMAGE ENHANCEMENT 
Once the beam sharpening techniques discussed above have been applied, image processing can 
further improve the output and generate more accurate and useful information. In order to 
implement these techniques, a high throughput processor is required such as the APG-70, APG-73 
or Hughes CIP. Several image processing techniques which have been widely implemented are 
presented below. 
1. Histogram Flattening 
In many images, detail in the darker regions is imperceptible. Signal processing is required to 
alter brightness or improve contrast. One method which works well is histogram modification. Both 
adaptive and non-adaptive histogram modification techniques exist. In principle, these techniques 
involve rescaling the image to some desired form. A common method of scaling is to find the 
average Value of the histogram and normalize quantized bands of pixels against it. By allowing fewer 
levels of quantization, higher contrast of features is obtained. This process, however, does result in 
an increase in quantization error. 
2. Edge Detection 
A common task in radar mapping is to determine the extent and position of extended targets 
such as airfields and roads. One method of performing this task is through the detection of 
luminescent edges or discontinuities between reasonably smooth regions. The only digital images 
which exhibit step edges are artificially generated test patterns and images. Due to low-pass filtering 
prior to digitization, any digital image resulting from optical or radar images of real scenes will have 
reduced edge slopes. For this reason, additional processing is required to identify and localize the 
boundaries of any extended objects. 
One of the primary methods of edge detection 1s differential detection. This can take the 
form of either first order derivative or second order derivative processing. In the first order case, 


some form of spatial first order differentiation is performed, with the resulting edge gradient being 
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compared to a threshold value. A gradient above threshold denotes the presence of a boundary. 
This process is demonstrated in Figure 20. [7] 

The simplest method of generating gradients is to compute the running difference of pixel 
magnitude along rows and columns. Diagonal gradients can be obtained by forming a running 
difference of diagonal pixel pairs. The running difference technique is highly susceptible to small 
fluctuations in object luminance. Object boundaries are also not well defined. Modifications do 
exist to better define edge locations. One such method which has provided good results is to simply 
take the pixel differences separated by a null. Luminance fluctuation problems have been limited 
through the use of weighted averages and moving windows. Larger windows improve the detection 
of edges in high noise environments, but the computational requirements increase rapidly with size. 
A sample C language program to perform first order derivative edge detection 1s presented in 
Embree [7]. This program uses convolution smoothing and vertical edge detection with a 3X3 
moving window called Sobel detection. 

In second order derivative processing, a change in the polarity of the second derivative 
signifies an edge. One implementation utilizes the Laplacian which is zero if the intensity is constant 
or changing linearly. A greater rate of change will result in a sign change at the point of inflection. 
The location of the zero crossing signifies the assigned location of the edge. A shortcoming of the 
Laplacian is that it produces a fairly noisy result. Extraneous edge points are often introduced, and 
especially strong responses result from corner reflectors. 

An alternate implementation Is to first estimate the edge direction. The second order 
derivative is then taken along the edge. The difficulty in this method hes in estimating the edge 
direction. First order detection techniques are sometimes used for this purpose. 

A second method of edge detection is model fitting. This involves fitting the observed data in 
a region of pixels to matrix mappings of possible edge types. Once a close fit is found, an edge is 


deemed to exist and the region of pixels is assigned the properties of the model selected. In this 
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Figure 20. 
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method, the match is made if the mean square error is below a threshold value. The model fitting 
method requires much more computation than the derivative detection methods. 

Pratt [6] provides an excellent comparison of these methods with several implementations of 
each. The implementations are evaluated with respect to three criteria. The first, good detection, 
refers to the ability of the system to detect edges without undue false detections. The selection of 
threshold mirrors the tradeoffs of probability of detection (P,) and probability of false alarm (P,) in 
target detection. The second, good localization, refers to the edge being designated as close to edge 
center as possible. The third, single response, states that only one peak above threshold should exist 
for each edge crossing. 

3. Line (boundary) Detection 

Similar techniques are used in the detection and refining of image lines. A line can be 
modeled as a series of closely spaced parallel edges. Weighted averages and model fitting for 
orientation are often used to limit susceptibility to noise. 

4. Thresholding and Centroiding 

This technique can detect and localize point scatterers as well as determine the extent of large 
objects. Thresholds are determined by local noise estimates. An m out of n detection scheme is 
employed to separate point reflectors from extended objects. If greater than m cells in an n cell 
window exceed the threshold, then an extended object is determined to exist in the space. If fewer 
than m threshold crossings exist, then a centroiding algorithm is applied at the peak crossing. This 
algorithm estimates the location and amplitude of the point object and designates it to the 
appropriate cell. The rest of the window is then replaced by the local background. An estimate of 
the location for a point scatterer is obtained by summing the weighted amplitudes of each cell and 


dividing by an unweighted sum. 


5. Median Filter 
Noise can exist in one of two forms. It can be contained in a separate portion of the 
frequency spectrum from the signal, or it can coexist with the signal across the spectrum of interest. 
In the first instance, linear filtering can be quite effective. In the second case, filtering to remove 
noise will also degrade the signal. An example of the second type of noise is "salt and pepper” noise. 
It often results from A/D converter problems or digital transmission errors. These bit errors result 
in impulses in an otherwise smooth sequence. Linear filtering to remove these “specks” in an image 
can result in blurred images and loss of high frequency information. An effective method of 
removing this noise is median filtering. 
Embree [7] describes median filtering as a simple four step process: 
1) A window of contiguous data is selected for each output point. This 
window can be reflected as adjacent time samples in a sequence, or as 
neighboring pixels in an image. 
2) The data is sorted by signal value, from high to low. 
3) The central value is selected as the median. 
4) This median value is used as the filter output. 
The advantage of this process is that sharp edges are preserved and not blurred as in averaging 
filters. A twist to this process is conditional median filtering. This technique avoids unnecessary 
corruption of the signal while maintaining the impulse removal capability of the previous filter. The 
process is as follows: 
1) Select window. 
2) Sort data from high to low. 
3) Determine median. 


4) The output of the filter is the median if the absolute value of the 
difference between the median and input is above a chosen threshold. 


Otherwise, the input is passed through as the output. 
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Median filtering is highly dependent on the type of sort used. In image processing, max or 
min values are often used to alter shrink or grow bright areas in an image. This is known as erosion 
or dilation. Selection of output can also be used to control contrast. 

Figures 21 and 22 [7] show the effect of median and conditional median filtering. These plots 
are the result of the program MEDIAN which is shown in Embree [7]. The "flat topping” evident in 
the figure is the result of median filtering a cosine wave. This is unwanted when filtering sinusoidal 
signals, but is acceptable for image row and column processing. It should be noted that the amount 
of “flat topping” is reduced through the use of conditional median filtering. 

Richards [8] demonstrates a further usage of median filtering to enhance noncoherent radar 
data. The algorithm diagrammed in Figure 23a [8] shows the combined usage of median and inverse 
filtering. A median filter is used to segment real air-to-ground radar data into two classes, potential 
point targets and background. The separation of the data overcomes the contradictory processing 
goals of targets and background, namely to enhance target spikes while smoothing background 
clutter. Once separated, the two classes of data are passed through deconvolution filters. This 
inverse filtering removes the known spatial antenna pattern from the data to yield a high resolution 
image. The resulting data is then recombined to yield an improved return. A lower order median 
filter is also used as the background component constraint in Figure 23a. [8] This was done to 
further control clutter spikes and ringing in the background which would otherwise be accentuated 
by the far lower diagrams. Figure 23b [8] shows the effect of this process on millimeter-wave data. 
Notice that the target response is sharpened in each case while the background has been noticeably 


smoothed. 
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Figure 21. Median filtering with sort lengths of 3 and 9 threshhold = 0. 
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Figure 22. Conditional median filtering with sort lengths 3, threshhold 
.l1, and sort length 9, threshhold .2. 
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Figure 23. An inverse filtering process utilizing median filtering (a) 
sampled data and filtered results shown below Cb 


IV. CONCLUSIONS AND RECOMMENDATIONS 

While the various methods of beam sharpening have been treated separately in this paper, in 
practice their implementation is often combined. This was demonstrated in the discussion of inverse 
filtering and the example from Richards [8] in which used median filters in the deconvolution 
process. In one Hughes system, the Hughes Advanced Discrimination Technique is used to detect 
the location of point scatterers, a modified inverse filter sharpens extended objects, median filtering 
is used to reduce impulse noise, edge detection algorithms are employed to identify roads and 
runways, and histogram flattening is implemented to improve image contrast and emphasize low 
contrast detail. When used in a system with high pulse compression, small mapping cells and 
improved resolution are obtained. 

Through the combined usage of these techniques, beam sharpening ratios of 6:1 can be achieved 
for isolated corner reflectors with a 6 dB signal-to-noise ratio (SNR). With a 20 dB SNR, a 25:1 
beam sharpening ratio can be achieved for corner reflectors and a 13:1 improvement for extended 
objects [5]. To illustrate the value of this capability, the following example will compare SAR and 
sharpened real beam resolutions for a fictional scenario. 

Assuming operation at 12 GHz, the border of X and K, bands, the system wavelength will equal 
.082 feet. If antenna length is assumed to be 10 feet, then the azimuth resolutions can be easily 
calculated. Stimson [1] gives the azimuth resolution of a SAR system as one half the physical 
antenna length. A comparison of resolution for real beam, sharpened real beam and SAR systems is 
shown in Table 6. From the unsharpened real beam results, it is not hard to see why these systems 
have been limited in application. The sharpened real beam columns, however, demonstrate "near 
SAR" quality. These examples do assume corner reflectors. The beam sharpening results for 


extended objects would only be approximately half as good. 
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TABLESG 
COMPARISON OF SAR AND REAL BEAM RESOLUTIONS 


RANGE SAR UNSHARPENED SHARPENED SHARPENED 
(in nmi) RESOLUTION REAL BEAM REAL BEAM (6 REAL BEAM 
dB) (20 dB) 





This level of improvement implies great potential for future mapping radars. Referring back to 
Table 1 on page 7, with a 20 dB SNR highways can be mapped to over 50 miles. "Road map" level 
of detail is available to over 20 miles. Taking into consideration that Scud launchers and other 
tactical targets are much larger than standard vehicles and are very well represented by corner 
reflectors, Table 6 demonstrates the capability to perform the mission proposed in the introduction 
of high-speed search and target localization. Further, these results were obtained by processing 
magnitude data only. In future efforts, these techniques may yield even better results through the 
processing of IQ data. 

Figures 24 through 29, provided by Hughes, demonstrate the effectiveness the techniques 
discussed in this paper. These figures provide a visual comparison of the results of these techniques 
and SAR imaging. This data was collected by an F-15 in tests on another project. Efforts and 
equipment dedicated to mapping in future tests may also yield further improved results. 

Resolution in range is equally impressive. As mentioned in Chapter II, bandwidth is the best 
determinant of performance. Since the 3dB bandwidth of a pulse is approximately the inverse of the 
pulse width, fine range resolution implies the need for high bandwidth. Continuing improvements in 
system bandwidth performance make it easier to achieve small cell size. 

Continuing efforts to increase signal bandwidth and improve azimuth resolution yield numerous 


opportunities for further research. Of particular interest in the future will be continuing efforts to 
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improve the extended object capability of mapping systems. Image processing can play a great role 
in this effort. With stand-off weapons and associated tactics evolving at an extraordinary rate, radar 


mapping will have a large role to play in any future conflicts. 
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beam sharpening. 





APPENDIX 
MATCHED FILTER OUTPUT AND THE AMBIGUITY DIAGRAM 
A. THE AMBIGUITY DIAGRAM 
One of the primary tools which will be used to evaluate these coding schemes is the ambiguity 
diagram. The ambiguity diagram represents the response of a matched filter to a waveform, and its 
doppler shifted versions, reflected from a point scatterer. 
1. Development 
The output of a matched filter can be represented as the cross correlation between a 


transmitted signal and the received signal. Neglecting noise, Skolnick [9] presents this as 
| S.(t)S (t-Tp')dt (A.1) 


where 
S(t) is the transmitted signal 
S(t) is the received signal. 
Tp’ is the estimate of the time delay 
The transmitted and received signals are assumed to be of the form 
S(t) = u(t)?" (A.2) 
S(t) = u(t-Tel? * Ot Oe) (A.3) 
where 
f is the carrier frequency 
f, is the doppler shift 


T is the time delay 
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Substituting these representations into (A.1) and simplifying by setting T=0, f=0 and - 
Tp'=Tp, the output of the matched filter becomes 
X(Tpf,)= | u(t)u (t+Tp)e?* (A.4) 
where 
u(t) is the complex modulation function 


|u(t)| is the envelope of the real signal 


With this representation, a positive f, denotes a closing target while a positive Tp denotes a 
target beyond the reference delay. The ambiguity function is the magnitude of this equation 
squared. 


The ambiguity diagram has a number of important properties as mentioned earlier. Several 


include: 
IX(Tpsfa)|7=1X(0,0) |?=(2E) maximum value (A.5) 
tee) cg he symmetry (A.6) 
i | |X(Tp fy) |?dT,df,=(2E) constant total volume (A.7) 


Eq. (A.7) is known as the waveform uncertainty principle and illustrates the problems and 

tradeoffs associated with waveform selection. 
2. Usage 

While the ambiguity diagram ts of limited use as a practical design tool, it 1s helpful in 
examining the limitations and utilities of particular classes of radar waveforms. By studying the 
ambiguity diagram for a particular waveform, judgements can be made regarding its suitability for 
various applications. 

Waveforms are selected to satisfy five major requirements. These are: 
1) detection, 2) measurement accuracy, 3) resolution, 4) ambiguity and 5) clutter rejection. All of 


these characteristics can be visually evaluated using the ambiguity diagram. The ambiguity diagram 
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allows the user to intuitively evaluate waveforms at a glance. General characteristics are presented 
clearly without intensive mathematical formulation, as demonstrated below. 


1) Detection is independent of the transmitted waveform and depends only on the level of 
energy which exceeds a chosen threshold at the origin. The maximum value of the ambiguity 
function occurs at the origin with a value of (2E)’. The threshold must be set below this value, 
but high enough to avoid any secondary peaks. 


2) Range accuracy is dependent on transmitted bandwidth (B) and doppler accuracy is 
dependent on the pulse width. The effects of these elements on the ambiguity function are 
shown in Figure 30. [9] Time delay (range) and frequency (doppler shift) are plotted on the 
horizontal axes with the ambiguity function magnitude plotted on the vertical axis. As will also 
be shown later, the volume under the ambiguity diagram remains constant. This is known as the 
radar waveform uncertainty law. Care must be taken in attempting to compress the diagram 
along either the range or frequency axis. (The waveform uncertainty principle forces a trade-off 
between range and doppler accuracy. Since the volume under the ambiguity function must 
remain constant, efforts to compress the diagram along one axis may result in a spreading of the 
diagram along the other. 


3) Resolution is closely related to measurement accuracy. A waveform with good resolution will 
always have good accuracy. A waveform with good accuracy, however, may not have good 
resolution due to the presence of sidelobes. Waveform resolution is displayed as two 
overlapping ambiguity diagrams displaced by the range and Doppler spacing of the targets. 
While waveforms with narrow spikes may provide excellent accuracy and resolution in the main 
peak, they may also have large sidelobes. Sidelobes from a large target or ground clutter may 
contain enough energy to overpower a smaller target's central peak. 


4) As alluded to above, peaks other than at the origin yield ambiguity in the measurement of 
range and doppler. Decisions with regard to pulse train length and PRF can be quickly 
evaluated by observing the presence and height of secondary peaks with respect to the main 
peak. 


5) Clutter rejection ability is visualized by insuring ambiguity function peaks do not exist in 
doppler or range bins with high clutter. 


3. Types 
There are five major classes of waveforms which give rise to three types of ambiguity 
diagrams. These types are the "knife edge", "bed of spikes" and "thumbtack". These are 
demonstrated in Figures 30. [9] 
The "knife edge", or "ridge", diagram is obtained from a single pulse of sine wave. Linear 


frequency modulation rotates this ridge through the fd,T plane. 
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Figure 30. The ambiguity diagram, ideal approximation and three major classes. 
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The "bed of spikes" results from a periodic train of pulses. The form of each of the spikes is 
dependent upon the waveform of the individual pulses. This is representative of pulse-Doppler and 
MTI systems. 

The "thumbtack" ambiguity diagram results from noise or pseudonoise codes. These codes 
include Pseudorandom, Barker and Polyphase codes. This thumbtack usually exists on a large 


pedestal of range-Doppler sidelobes which can be quite significant. 
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B. AMBIGUITY DIAGRAM PLOTTING PROGRAM 


pe me eee ie ee 
— <_< ee eee eee ee es ees ee eee es ee eee ee ee ee es es es es es es es ese ei ee eee 


ambigUus Michael D. Anderson 


= A menu driven plotting program to draw ambiguity = 
= diagrams for various pulse compression code families. = 
= AmbigUus is the main, calling menu subroutines and = 
= performing the mesh plot. = 


clear 

elg 

CODE — lL, 

alte ae (Calls subroutine to provide instructions) 
waveform (Calls subroutine to route user through menus. 


Menu programs implement selected codes and 
return for phase incrementation and plotting) 


oP oP olf ol? 


o 


* (Scale the ambiguity diagram by selecting max reach of the 
doppler axis.) 


S=input ('Enter max doppler phase shift across the pulse in 
radians: ') 


l=length (CODE) ; 


EOi =e or % (Implement the iterative phase 
step=(s*2)/(50-1); % shift across the code for each 
shift=(-s)+(k-1)*step; % increment along the doppler 

6 axis.) 
jae yeails I< 
v(j) =shift*j+CODE(j); 
end 
Ua=1* CODE; 
Va=i*v; 


ub=exp (ua) ; 
vb=exp (va) ; 
QUE DUE —XeGien ule, vom 
x=abs (output) ; 
X=x.*2; 
Ae. =x 

end 

mesh (A) 


(Digitalizes the phase code values) 
(Digitalizes the shifted phase code) 
(Performs the autocorrelation) 

(The matched filter output) 

(The ambiguity function is the 
Square of the filter function) 


oe o\e o\e o\e ol? ole 


ole 


(Plots the 3-D ambiguity diagram) 
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omc mcr rm wc mc mw wm ww mm mw mc ww cr we eee i 
NN _ _ ____.__.__._.._........._...._...__..._...._..._ ee  ___._._  ___ _____  ________ 


= intro Michael D. Anderson 


= Greets users and instructs them of the purpose = 
= and usage of the program. = 


ee mc wm cc wm mw cr cw ww wc cw mw cw ww ee i 
ec cc cr cc cr cr cc cr cc cr ec ce cc ce cc ei 


Welcome to the program ambigUus. This program will enable 
you to view the ambiguity diagram for a number of your favorite 
waveforms. 

The ambiguity diagram is a 3-D plot showing the response of 
a matched filter receiver to reflections from point scatterers. 
Range and frequency shifts are read along the horizontal axes, 
as the matched filter response is plotted along the vertical 
axis. 

This plot can be used to evaluate the critical properties 
of detection, range and doppler accuracy, range and doppler 
resolution, range and doppler ambiguity, and clutter rejection. 


(press any key to continue) 
echo off 


pause; 
end 
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(qm pte ccc ccc ccm mcr cm crm mmc ec ccc cc ee ee ee ee ee eee ee 
mmc ccc cc cmc cmc rece mcm cc ec cee ce ccc ce ee ee er ee ee ee eee ee eee ee eee ee ee ee ee 


= waveform Michael D. Anderson - 


= A menu program allowing the user to select which = 
waveform to implement. = 


em mcr mcm mm mmm cr mmm cc rm mm mc mr eee ee eee ee 
mete emer cc cm ccc ccc ccc ccc cc ccc ce cc cc cc cc cc ee ee ee ee ee ee ee ee ee 


echo on 
Two major classifications of radar waveforms are pseudo- 


random binary codes and polyphase codes. Both can be modeled 
with this program. 

Pseudo-random binary codes consist of Barker, Compound 
Barker and Complementary codes. 

Polyphase codes consist of Frank codes and the P-series 
codes. These are digital representations of LFM and step-chirp 
waveforms. 


Barker code 

Compound Barker code 
Pseudorandom code 
Polyphase code 

Frank code 

P1 code 

P2 code 

P3 code 

P4 code 

quit program 


OWoOWMO NAM PWNEH 


oP of ofP off ofP o'P aff OP OP O'P OfP O'P O'P O'F O1P O19 O'P OP OLP O'P OP OP O'P O'P O\P OP O'P OP 


echo off 
n=input ('Select the waveform type to process: ‘'); 


1f n==1 
barkform 
elseif n== 
compound 
elseif n== 
randform 
elseif n==4 
polyform 
elseif n==5 
frnkform 
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elseif n==6 


piform 
elseif n==7 
e25Orm 
elseif n== 
po onrm 
elseif n== 
p4form 
else 
death 
end 
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‘mmm ems mmm cc meee ccc ce ecm cr crc ccc tm ce wt tt cm cr cr cr cr wr crt ww wt ee ee ee ee 
em ccc eee ccc cmc cmc crc crete tte ttc cr crc cc cc cc crt cr crete tee ee eee eee eee ee ee eee ee ee 


= Josie izieonant Michael D. Anderson = 


the vector CODE with the code phases. 


emcee mmc cmt ttre ttre tte ttt cc cr ct tert tet tte te cr cr teste tte tc tet ee ee eee ee ee 
| Sas my | cmc cme cece eee eee ccc crc cr ctr emt ete tem te ce te ete ete te ee ee ee ee eee ee 


6 = Allows user to select code length to plot and loads = 


The Barker codes are a special set of binary phase codes 
which yield equal sidelobes after passage through a matched 
filter. These codes have pulse compression ratios equal to 
their length. The (+) and (-) in the sequences represent 
phase shifts of zero and pi radians. The known Barker codes 


are as follows: 


selection code # code sequence 
. 2A foe 
é 2B fant 
3 3 ++- 
4 4A ++-4 
5 4B +++- 
6 5 $++~-4+ 
Zi 7 Ng 
8 ca +++---+--4- 
S) is t+t+¢4+--4¢4-4-4 
SCO Off 


n=input ('Enter your choice, selection # '); 
if n==1, 

GObE— Cr oi, 
elseif n==2, 

GCOBE=[0 0])- 
elseif n==3, 

CODE Os OF pall. 
elseif n==4, 

CODE=[0 0 pi 0]; 
elseif n==5, 

COD i=" |10) 18) 10 jes. |= 
elseif n==6, 

CODE=[0 0 0 pi 0]; 
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elseiff n==7, 
Sout 0 OU Umm. O pi) - 
elseif n==8, 
SOPE— 10 100 pis pty p2 0 papi 0 pil; 
elseif n==9, 
COUE— | O semOmOme soll On Cup 1 CaoieO)| | 
end 
end 
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mw mc mc mw mm ww mw mw mc ww cm cm rw eee 
— me ae ee eee eee ee see oe cee eee ce ce ce ce es ee es es es a i es es es ss es es es es es es ei ees es es es es es es es ei ee ee eee ee 


= compound Michael D. Anderson 


= compound barker codes. Originated from "combined.m" 
= with Capt. Paul Ohrt,) RGA@ ter Fer49 707 oe 


emcee cm ww rr ws we we ee i 
8, ee —_ |_| 


2 
°o 
2 
°o 
o 
o 
oe Allows user to select inner and outer codes to generate = 
6 
2 
o 
o 


clear 


BARKERG— (0, pal 

BARKER3=[0,0,pil; 

BARKER4=[0,0,pi,0]; 

BARKER5=[0,0,0,pi,0]; 
BARKER? = [Uy 090 y oan Otero ann 

BARKER? 1=([07; 0; 0 709%, 09 ou) 0 pi pte eon 
BARKER13=[0,0,0,0,0,pi,pi,0,0,pi,0,pi,0]; 


oe 
echo on 
This subroutine generates Compound Barker codes. 
The known Barker codes have lengths of n= 2,3,4,5,7,11 or 13 
You need to choose the length of the inner Barker code, "m" 


and outer Barker code, "n". The length will be m*n and the 


PSL will be at the level of the shortest Barker code used. 


OP o\9 a\9 ofP ofP oP aff o'P O'P OP 


echo off 
m=input ('ENTER the # of the inner Barker Code to use, m= ') 
n=input ('ENTER the # of the outer Barker Code to use, n= ') 


=e ™ Oe 


° 


% Set selected inner Barker Code 


if m== 

in _CODE=BARKER2; 
elseif m== 

in _CODE=BARKER3; 
elseif m==4 

in _CODE=BARKER4; 
elseif m==5 

in_CODE=BARKERS; 
elseif m== 

in _CODE=BARKER7; 
elseif m==11 

in_CODE=BARKER11; 
elseif m==13 

in_CODE=BARKER13; 
else 


RZ 


° 


compound 6 Starts sub program over if improper "m" entered 
end 


% Set selected outer Barker Code 


if n==2 

out _CODE=BARKER2; 
elseif n== 

out _CODE=BARKER3; 
elseif n== 

out _CODE=BARKER4 ; 
elseif n== 

out_CODE=BARKERS ; 
elseif n== 

out _CODE=BARKER7; 
elseif n==11 

out CODE=BARKER11; 
elseif n==13 

out _CODE=BARKER13; 


else 

compound 6 starts sub program over if improper "n" entered 
end 
6 Take Cosine of phase changes and develop "m in n" BARKER CODE 


in cos=cos (in CODE) ; 
out cos=cos (out_CODE) ; 


CODE={]; 

fPomecount=1:n 
imIe—out cos (count) .*in cos; 
CODE=[CODE, INC] ; 

end 


CODE=acos (CODE) ; 
end 
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= randform Michael D. Anderson = 
= Allows user to select code length to plot fora = 
= polyphase code using shift register implementation. = 


= Program originated from "random.m" developed with = 
= Capt] Paul Ohnt RCA for HC4970rw 1992 = 


When a large pulse compression ratio is required, 
pseudorandom codes are often used. These approximate noise 
modulated signals with "thumbtack" ambiguity functions and 
excellent range and doppler resolution. 

A commmon method of generating these codes is the modulo 2 
adder or exclusive-OR gate shift register decribed in Skolnik 
ch 11. That method is used here. The only restrictions are 
that the initial conditions of the register can not all be zero 
(or the output will simply be a sequence of zeros), and that (p) 
can not equal (q) (both shift inputs can not be taken from the 
Same bit). The length of the resulting code will be —- 

(hit any key to continue) 
pause; 
Che 
n=input ('ENTER the length of the shift register (# of bits): '); 
p=input ('ENTER the first feedback bit: # '); 
q=input ('ENTER the second feedback bit: # '); 


echo on 


To enter your initial conditions, simply enter a string of 


ones and zeros. 


oP oP olP olf olP ol 
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ex. Shift register length = 7. Initial conditions = 1001100 


of off ol? 


eecno oft 


flag=1; 

while flag=1 
me maput('Enter your string of length [n]: '); 
tielength (IC) =n, 


flag=0 
end 
end 
SODE=[IC]; So ANINtial Stare pont 
MLS=2“n-1; % maximal length sequence 


Pemeeount = (n+1):MLS 


XOR=CODE (p) | CODE (q) ; % OR function from Matlab 
if CODE(p) &CODE(q) ==1; 

XOR=0; %* to ensure XOR conditions met 
else 

POR =x ORK? 
end 
CODE= [XOR, CODE] : 6 update code values from shift 


end 
CODE=acos (CODE) ; 


end 


cD 


—— eee ccc ccc cc cc mm cr crm cm cc cw ww wm eee eee ee 
— er yee cc ccc ccc ec emer cr creme cr ce ese ec ce cee eee ee eee ee ee eee ee eee ee 


= pels tem Michael D. Anderson = 


= polyphase sequence and loads the CODE vector. Program ~ 
= originated from polyphas.m with Capt. Paul Ohrt RCA = 
= “form BC 4970732. = 


— crc crm crc cr cr cr cr cr cr cr cr cc cr cr cr cr wr cr mc crm cm ee eee ee ee 
—— cm ccc cece ccc ccc ccc cc cre ec cmc crc crc crc cc cc ee ee ee ee ee see ee ee ee eee eee eee eee ae ae a oe 


a Allows user to select code length to plot for a simple = 


6 This subroutine implements a POLYPHASE coding 

6 sequence. The POLYPHASE code developed here 1S a sequence 
6 between 0 and pi inn increments. 

6 You must specify how many increments are to be used. 


n=input ('ENTER the # of increments the Code is to use, n = '); 
6 (Generate polyphase vector) 

k=0: (n-1) ; * (increment vector) 

phi=pi.*k.*2./n; 

CODE=cos (phi) +1.*Sin(phi) ; 


end 
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ti i ______  —  _ _ e 


= miictorm Michael D. Anderson = 


= the vector CODE with the code phases. From "frank.m" = 
= generated in 1992 for EC 4970. = 


o.= Allows user to select length of a frank code and loads 


clear 

eile 

echo on 

8 This subroutine implements Frank Coding. It is equal to 
: Bamgker coding in PSL, but 2s not Timited in length. 

% Frank coding is a step-wise approximation of a 

: linear FM or "chirp" code. To use this coding scheme, 

: enter the number of phase subdivisions N. The PCR 

z will be N squared. 


Nemmewe('’N = '); 
step = 2*pi/N; 


Pome = 1:N 
nee | = 1:N 


fount = (k=-1)*N + 5; %* Generate the Frank Code Matrix 
eeeouttt) = (k-1)* (7-1); % Store in vector form 
end 
end 
stepx = step*x; %$ combine the Frank Code Matrix with the actual 
$ phase increment to generate the phase sampling 


CODE=cos (stepx) +}.*Sin(stepx) ; 


end 
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i LL ee 


= pit Gig Michael D. Anderson = 


= the vector CODE with the code phases. Derived from = 
= Ernktormmem 1992) ~ 


em ccc cr mc cmc ee me cm cr cc ccc cc ccc crc mm ccm cc ec ee ee eee ee eee ee eee eee ee eee oe 


6 = Allows user to select code length for Pl code and loads = 


echo on 
This subroutine implements Pl coding, a Frank code 


rearranged to avoid bandwidth limitations. The ambiguity 
diagram will be almost identical to the Frank code since both 
are derived from the same waveform. 


Enter the number of phase subdivisions N. 


OP olf o\P off olf off olf oP O'P OP? 


echo off 


N=input('N = '); 
step = pi/N; 


fei — aN 
Ome aaa sea 


count = (k-1)*N + j; 
x(count) = ((j-1)*N+(k-1)); 
x(count) = x(count) * (N- (2*j-1)); 
end 
end 
stepx = (-1)*step*x; 


CODE=cos (stepx) +].*sin(stepx) ; 


end 
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- p2form Michael D. Anderson 


= Allows user to select code length for P2 and loads 
the vector CODE with the code phases. 


echo on 
This subroutine implements P2 coding. P2 coding 


1s a real vice complex version of Pl. Results should 
be indistinguishable from Frank and Pl codes. 


Enter the number of phase subdivisions N. 


OP o'P al? ol? Ol? ol? O19 oP 


echo off 


feomeut('N = '); 
step = (pi/2)*((N-1)/N) ; 


femek = 1:N 
mon | = 1:N 


count = (k-1)*N + Jj; 

x(count) = step- (pi/N) *(k-1); 

x(count) = x(count) * (N+1-2*74); 
end 


end 
CODE=cos (x) +} .*Sin(x); 


end 
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emer wr ww ww ww www tw ee ie 
mercer crc mre tere tec ee eee eee ee ie 


emmy rr rrr crc rw ww i 
‘cece rrr crcrrcrrecrcmcrrcrc ee e eeeeee ee 


Michael D. Anderson 


= the vector CODE with the code phases. 


ecm ee ee ee eee ee ee 
ce cee ee ee ee ee ee eee ee es es es es eo ee ee ee es ee es es ees es es es ee es es es es es ee es ee ee ce ee eee ee ee ee 


echo on 


oP of© olf olf off of olf olf off of? 


echo Gif 


N=input('N = '); 
Step = pi/N; 


OG 


LJ 
ll 
er 


Ny 
jake 


stepx = step*x; 
CODE=cos (stepx) +].*Sin(Stepx) ; 


end 
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chirp as the Frank, Pl and P2 codes. 


frequency modulated (LFM) signal. 


6 = Allows user to select code length for P3 and loads 


This subroutine implements P3 coding. This code 
is derived from a linear chirp rather than the step 
Because of this, 


it will show the better doppler tolerance of the linear 


ccm mmm ccm wm cw ww ie i 
a a a ert i i i ei 


= p4form Michael D. Anderson a 


= Allows user to select code length of a P4 and loads = 
the vector CODE with the code phases. = 


emcee wr cr rc cc cw ec ww we cw ww ww ww ww ie i 
cep cc ccc cr ccc cc cme ce ec cc ise 


This subroutine implements P4 coding. P4 codes 
are simply P3 codes rearranged to insure that the 


greatest phase shift between code bits occurs at the 


identicle in conncept to the reworking of Frank codes to 
Pl to remove the bandwidth limitations. 


© 

6 

2 

° 

2 

fo] 

a 

° 

2 

o 

% edges of the sequence rather than at the center. This 
6 

2 

° 

2 

° 

2 

° 

2 

°o 

% Enter N, the number of phase subdivisions to use. 

2 
°o 


N=input('N = '); 
Seep. —eoi/N; 

moog = 1:N 
pe () ) 
mt) 
end 


| 
-_— 
Ke 
| 
be 
~~ 
~e 


CODE=cos (x) +j].*sin(x) ; 


end 
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